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1. Introduction 

The successful LHC run during 2010 - 2012 has resulted in the accumulation of large samples 
of data of the Standard Model processes like the heavy gauge boson production and the decay 
to lepton pairs along with the announcement [1] of the Brout- Englert-Higgs (BEH) [2] boson 
candidate. The era of prediction of QCD processes at sub-1% precision tag is upon us. In order 
to obtain this desired level of accuracy the infrared (IR) improved DGLAP-CS [3, 4] theory [5, 6] 
realization with HERWIRI1.031 [7] was done by implementing the set of IR improved DGLAP 
kernels in HERWIG6.5 [8, 9]. It has been argued that this process allows better than 1% theoretical 
precision [10, 11]. The comparisons previously done with ATLAS [12] and CMS [13] data are 
encouraging. In the following we extend this to the LHCb [14] results. 

Any precision theory should cover the entire observable phase space for hard processes so that 
it can completely exploit the entire LHC data. Eor the single Z/y* production and the subsequent 
decay to lepton pairs, the LHCb probes the regime of pseudo rapidity given by 2.0 < T] < 4.5 
which is different from the regimes probed in the data in Refs. [5, 6]. Thus the LHCb provides an 
oppportunity for a different check. It is to be noted that the cuts on the lepton transverse momenta 
were similar in all the data. 

In what follows we will discuss the comparison of the IR-improved and unimproved parton 
shower(PS) MC predictions, with the MC@NLO [15] PS/matrix element (ME) matched exact 
t^{as) correction, with the LHCb data on the Z/y* rapidity, pr and 0* distributions. The variable 
0* was introduced in Refs. [16] in order to overcome the difficulty of measuring pj spectra for the 
lower regime of pj. The definition of the variable will be provided in the sections to follow. 

In the next section we will give a brief review of the theory of Exact Amplitude-Based Re¬ 
summation. In Section. 3 we show the comparison with the LHCb data. We will then discuss the 
attendant theoretical implications. 

2. Precision QCD for LHC 

We start the discussion in this section with the fully differential representation of a hard LHC 
scattering process in order to link the experimental results with the theoretical predictions: 



( 2 . 1 ) 


where {F, } and ddres are the parton densities and the reduced hard differential cross section respec¬ 
tively, the later of which has been resummed over all large EW and QCD higher order corrections. 

Eor both the resummation of the reduced cross section and the evolution of the parton densities, 
the defining formula is identified as 



Pi ^2 


( 2 . 2 ) 
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where dOres is either the reduced differential cross section dOres or tho evolution rate associated 
to a DGLAP-CS [3, 4] kernel involved in the evolution of the and where the new YFS-style 


[17, 18] non Abelian residulas Pn,m{ki ^have n hard gluons and m hard photons. 


Here the final state has been shown with two hard partons with momenta p 2 , qi specified for a 
generic 2/ final state. The infrared functions SUMjr{QCED) and Dqced are defined as 


SUMmiQCED) = 



(2.3) 


where K^ax is a dummy parameter and nothing depends on it. The following have been introduced 
in (2.3) 



anls 

^QED' 


(2.4) 


^QCED = ^QCD^^QED' 


The superscripts nls here denote that the infrared functions are DGLAP-CS synthesized. We stress 
here that in the formulation of equation (2.2) the entire soft gluon phase space is included. The new 
non-Abelian residuals j8 allow rigorous PS/ME matching via their shower substracted analogs. In 
equation (2.2) we make the replacement 



where all effects in the shower associated to the {F,} are removed from The MC@NLO 
differential cross section can be written as 


dOMC@NLO^ [B-\-V + j {RMc — C)d^E\d^B[^Mc{^)A- J {RMc/B)AMc{kTd^R)] 

+ (F — RMc)Ac{kT)d^Bd^R 

where the Sudakov form factor, which as usual represents the no emission probablity, is 


(2.5) 


B 


AMc{PT)=e^ 


0(kT(^B,^R)-PT)] 


( 2 . 6 ) 


B is the Born distribution, V is the regularized virtual correction, C is the corresponding counter 
term required in at exact NLO, R is the respective exact real emission distribution for exact NLO 
and Rmc is the parton shower real emission distribution. We show in Ref. [8] how it is realized via 
eqn. (2.2). 

The resulting new resummed kernels, yield a new resummed scheme for the PDF’s and 
the reduced cross section: 




(2.7) 


Pg,{z)^If/{z)^CFFYFs{Yq)e 
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The new scheme gives the same value for a with improved Monte Carlo simulation. Here the YFS 
infrared factor is Fyfs = +a) where Ce is the Euler’s constant. 

The new scheme has improved MC stability. In Herwiri 1.031 there is no need for an IR cut-off 
parameter ‘kf. The degrees of freedom below IR cut-offs that are dropped in the usual showers are 
included in the Herwiri 1.031 showers. We note that the difference in the showers starts in 
the new kernels agree with the old kernels at ^(a^). 

3. Consistency Checks 

MC HERWIRI1.031[7] is the first realization of the new IR-improved kernels in Herwig6.5 
[19] environment. In Refs. [7, 8] it was copared with the data of ATLAS and CMS on the single 
Z/t* productions and subsequent decay to lepton pairs. 

The MC@NLO/HERWIG6.510 simulations require a PTRMS= 2.2GeV to give good fits to 
both sets of data whereas the MC@NLO/ HERWIRI 1.031 simulations give good fits to both data 
sets without such an ad hocly large PTRMS [7,8]. 

We now move on to do the comparison for the LHCb data. Eig.l shows for the single Zj'f 
production at the LHC the comparison between the LHCb rapidity data for the e^e~ channel and 
the MC theory predictions. Fig. 2 shows the comparison with the LHCb rapidity data for the li~ 


LHCb Vector boson rapidity | 



Figure 1: Comparison with LHCb data: LHCb rapidity data on (Z/7*)production to e^e pairs, 
the circular dots are the data, the green(blue) squares are MC@NLO/HERWIG6.510(PTRMS = 2.2 
GeV/c)(MC@NLO/HERWIRI1.031). The green triangles are MC@NLO/HERWIG6.510. These are other¬ 
wise untuned theoretical results. 

channel. These results should be considered from the perspective of our analysis with the data from 
ATLAS and CMS in Refs. [7, 8]. What we found is that the IR improvement in HERWIRIL031 
allowed it to give a better jd.of to the ATLAS and CMS data. It did not need a large intrinsic 
value of the PTRMS. From already existing models of the proton [20] PTRMS ^ 0.4G^y/c, which 
is also what is indicated by the precociousness of Bjorken scaling [21, 22]. The unimproved results 
from HERWIG6.5 required PTRMS = 2.2GeV/do get similar /d.o.f values for the pr spectra. 
However, for the rapidity data such a large value of PTRMS was not required. The LHCb data 
comparison in Figs. 1 and 2 show a similar result. The values of x^/d.o.f are 0.746, 0.814, 0.836 
for the respective predictions from MC @NLO/HERWIRI1.031, MC @NLO/HERWIG6.5 (PTRMS 
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Figure 2: Comparison with LHCb data: LHCb rapidity data on (Z/y*) production to bare p~ pairs, 
the circular dots are the data, the green(blue) squares are MC@NLO/HERWIG6.510(PTRMS = 2.2 
GeV/c)(MC@NLO/HERWIRIL031). The green triangles are MC@NLO/HERWIG6.510. These are other¬ 
wise untuned theoretical results. 

= 0) and MC@NLO/HERWIG6.5 (PTRMS = 2.2 GeV/c) for the data and are 0.773, 0.555, 
0.537 for the respective predictions for the data. The values are acceptable for all three 

calculations. 

When we turn to the transverse momentum degrees of freedom, the situation is different from 
the Refs. [7, 8]. We start with the 0* data of LHCb [14]. The definition of this new variable is 

(j>* = iw{(l)acop/'2.)\J 1 -tanh2(AT]/2) (3.1) 

where, At] = 77 “ — T]+ where t]“ and 17 + are the respective negetively and positively charges lepton 
pseudo rapidities and ^acop = ^ — A0 with A0 = 0i — 02 which is the azimuthal angle between the 
two leptons. It is to be noted here that this variable is not the same as pj but it is correlated with it. 

Fig. 3 shows only the MC@NLO/A results, for A= HERWIG6.5 (PTRMS=0), HERWIG6.5 
(PTRMS =2.2GeV/c) and HERWIRI1.031(PTRMS=0) in comparison with the data for the LHCb 
0*. We find the respective jd.o.f in this case are 1.2,0.23,0.35 for the MC@NLO/HERWIRI 1.031, 


I LHCb Vector boson rapidity | 
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Figure 3: Comparison with LHCb data on 0* for the p^p channel in the single Z/y* production at LHC. 
The notation is the same as Fig. 2 
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MC@NLO/HERWIG6.5 (PTRMS = 0), MC@NLO/HERWIG6.5 (PTRMS = 2.2 GeV/c) simula¬ 
tions. We see that all three calculations give an acceptable fit but the MC@NLO/ HERWIG6.5 
(PTRMS=0) gives a mildly better fit than MC@NLO/ HERWIG6.5(PTRMS=2.2 GeV/c). 

In order to be consistent with the comparisons with the CMS and the ATLAS data we turn 
to compare the pr spectrum of the LHCb data with the MC theory predictions. Eig. 4 shows the 
corresponding comparison. 


[LHCb Vector boson | 



Figure 4: Comparison with LHCb data on pj for the p^p channel in the single Z/y* production at LHC. 
The notation is the same as Fig. 2 

The plot shows the comparison between the three MC@NLO/A predictions and the LHCb 
Pt data, where A = HERWIG6.5(PTRMS = 0), HERWIG6.5(PTRMS = 2.2 GeV/c) and HER- 
WIRI1.031, and as usual we always set PTRMS = 0 in HERWIRI1.031 simulations. The respective 
/d.o.f are 0.183, 0.103, 0.789 respectively. We see that all three calculations give an acceptable 
fit to the data, with a very mild indication that HERWIG6.5(PTRMS=2.2 GeV/c) gives a better fit 
than HERWIG6.5(PTRMS=0 GeV/c). We conclude that, when looking at the data on single Z/']Z 
production at CMS, ATLAS and LHCb , HERWIRI1.031 gives a good fit to the analyzed data 
without the necessity of ad hocly large intrinsic PTRMS. One of us (A.M.) thanks Dr. Kenneth T. 
Wilkins, for the kind hospitality of the Baylor College of Arts & Sciences. 
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